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Abstract 

The present study is both a reanalysis and an extension of the approach initiated by Eremin and 

o 

CN ' Voloshin(referenced in the text). We attempt to interpret here the rapidity-spectra of the various 

^ I 

■ particles produced in both Pb + Pb and Au + Au colhsions at CERN-SPS and RHIC-BNL. The 

O 

study made here is wider in scope and is more species-specific than the earher ones for which the 
results obtained here have been compared with those suggested by some previous works based on 
HIJING, VENUS etc, at various centralities. The study reconfirms that the constituent parton 
picture of the particles provides a better and more unified description of the rapidity-density yields 
for the various secondaries, including some light cluster particles like deuteron even in heavy ion 
collisions. 
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I. INTRODUCTION 

In the recent past Eremin and VoloshinQ] proposed that both nucleons and individ- 
ual nuclei could be considered as superposition of the constituents (of the particles) called 
partons/quarks/valons etc. Normally it is assumed that the nucleons are built of three 
such constituent partons and the mesons are composed of two of them. Such a constituent 
picture helps us to interpret the deep inelastic collision processes involving lepton-hadron, 
hadron-hadron and hadron-nucleus collisions. Much earlier, Bialas et al studied to explore 

n 

the hypothesis, within the framework of the additive quark model(AQM)[2], that the con- 
stituent partons(quarks) provide the universal elements not only of nucleon-nucleon and 
nucleon-nucleus reactions but also nucleus- nucleus interactions at high energies; this study 
ended up with surely an affirmative indication to such a possibility. But, as the focus was 
on different issues,the central theme of this work deviated from what it is intended to be 
here. We will simply look into the behaviour of some chosen observables studied here for 
nucleus-nucleus collisions alone from the viewpoint of such partonic constituent pictures. 
The centrality- dependence of multiplicity-density offers a very fundamental observable and 
we present here our studies on production of pions, kaons, protons and antiprotons sepa- 
rately in both Pb + P6(including deuteron production for this collision) at CERN-SPS and 
Au + Au collisions at RHIC-BNL at ^/s^n = 130 and 200 GeV respectively. Besides, very 
recently the data on deuteron production in Pb + Pb collision have been obtained, for which 
they have also been taken here into account. This apart, we will also compare the perfor- 
mance of the present approach with those obtained on the basis of some standard- version 
programmings like HIJING, VENUS etc which are essentially much more complicated than 
the one dealt with here. Our objective here is to examine the efficacy of this approach in ex- 
plaining some relevant data on production of these very important kind of secondaries in the 
few now available high energy nuclear collisions. The approach was also successfully tested 
for production of neutral particles, like photons, by Netrakanti and Mohantyls^ in the very 
recent past and was found to be in accord with the contention of Eremin and Voloshin^. 
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II. THE OUTLINE OF THE CALCULATION OF THE NUMBER OF PARTICI- 
PANTS 



In course of our calculations for the number of participant nucleons, denoted by Nn-part, 
and the number of participant partons, denoted by Ng_part we follow exactly the ways 
adopted by Eremin and Voloshin|]| and Netrakanti and Mohantyjsl, wherein a Woods-Saxon 
nuclear density profile is taken into account. This is given by 

''^^''^ = l + exp[{r-R)/d] 
where Uq = 0.17 /m'^, R = {l.UA^/'^ - 0.86^-1/^) ^ = o.54 fm. 

The number of participant nucleons (A^„_part) for a nucleus-nucleus ( A + i?) collision at an 
impact parameter b is calculated using the expression [l|, 3|, 

Nn-part\AB = J T^{s) {l - [ 1 - ^l^M^ ]B} 

+ J £sTnis-b){l -[1 - ^-^^]^}, 
where A and B are the mass numbers of the two colliding nuclei, T{b) = 
dznAiVb'^ + -2^) is the thickness function, and a]^^ is the inelastic nucleon-nucleon cross 
section. The number of participant partons is also calculated in a similar manner by taking 
into account the following changes of the facts related to physical realities, viz., (i) the den- 
sity is changed to three times that of nucleon density with = Suq = 0.51/m~^; (ii) the 
cross sections are changed to cTgq = cr]^^/9 and (iii) the mass numbers of the colliding nuclei 
are changed to three times of their values keeping the size of the nuclei same as in the cases 
of iV^^^lllQ. 

In the same way, the number of participant partons in a PP or PP collision can also be 
calculated by using A = 3 and 5 = 3, and considering nucleons as hard spheres of uniform 
radii 0.8 /mj^. 

For the present work we use cr]^^ = 30 mb in the range of cm. energies ^/snn = 4-17.3 
GeV as suggested in Ref.jSl. For higher energies, we first fit the data on total cross-sections 
for PP and PP collisions y] at energies beyond ^/Jmn = 30 GeV with the expression Q] 

= (^^•^^o.3)[i+ ^'-^ ^^^-^ ]+(i.i^o.i) m (^)+(o.i9^o.oi) 

(3) 



3 



with /ndf = 171/47, where ndf = No. of data — No. of Parameters. We also fit the 
elastic cross-section data for the same interactions with a similar type of expression which 
is given by, 



- (6.2±0.2,|lH-iH±i|pi,_(0.4T±0.04)ln(^)^(0.n±0.01),n^,^: 

(4) 

with /ndf = 53.14/18. The fits along with the data are shown in Fig.l. Hence, one can 
obtain the energy-dependence of inelastic cross-section(o"]v^^) for nucleon-nucleon interaction 
by subtracting eqn.(4) from eqn.(3). The values of cr^Xf) obtained in this manner, at some 
energies are given in Table-|l] The obtained values fall within 6% error of those suggested 
in Ref.^j. But a point is to be noted. Though some sort of justification for proposing the 
nature of the total cross section in the form of expression (3) given in Ref.j7|, none such 
explanation for expression (4) is possible. This is simply assumed with a view to providing 
only a best fit to the available data on elastic cross section shown in Fig.l and extracting a 
workable phenomenological form of expression for the nature of inelastic cross section, cr]^^. 
In fact, these twin relations help us to arrive at the usable values of inelastic cross sections 
needed for our necessary calculations of the relevant observables in a systematic manner 
even at energies for which no measured data on cross section values are available. 

III. RESULTS AND DISCUSSIONS 

The values of the number of participant nucleons and those of the participant partons for 
three nucleus-nucleus collisions at various energies and different centralities are presented in 
the Tables (jlTl - II Vj) . A comparison of the results obtained on the basis of the present work 
with those of NA49|8j, WAQsQ and PHENIX 1^, 
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12l | groups is made in Fig.-2(a) to 



T?- -4- f 1 <K-pa%f'^°''''> j.rNA49/WA98/PHENIX t, • . A 

Fig--2(d), m terms of 1 - na49/wa98/phenix ^ versus N^.pJt • K is observed 

that the agreements in cases of Au+Au collisions are modestly fair, whereas for Pb+Pb 
interaction the disagreement is quite strong and prominent. This discrepancy could, for the 
present, be attributed to the much less cm energy in lead-lead reaction which falls roughly 
at 17.2 GeV only. 

The integrated yields of the identified hadrons at central rapidity region produced in 
Pb + Pb and Au + Au collisions at different centralities are presented diagrammatically in 



Fig.3-Fig.5. In subfigures labelled as (a) and (b) of each category, the experimental data 
on integrated yields are normalized by half of the number of participant nucleons(A^„_part) 
while those in subfigures labelled as (c) and (d) are normalized by half of the number of 
participant paitons{Nq_part) ■ As we are here to make comparison between these two cases, 
we should use values of both the observables obtained from a single model. And that is why 



we put into use the values of < Nn-part > obtained on the 



IG 



3asis of the present work,instead 
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of the values indicated by NA49[8[ or PHENIX groups 

Now, to make a comparison between these two cases we fit the data in Fig.3-Fig.5 by two 
phenomenological expressions which are given by. 



1 dN, 



ch 



0.5 < 

—part 

and 



dNr. 



a A^^-paH (5) 



0.5 < N,.part > dy ^ 
where a, b, a and (3 are four constants. The fitted values of these parameters are given in 
Tables (fVTl Villi and depicted by solid curves in the figures. As could be seen from Fig.2 that 
the most peripheral values of Nn-part for Pb + Pb collision at Eiab = 158^4 GeV and Au + Au 
collision at y^J/vjv = 200 GeV obtained by the present model show larger discrepancies with 
respect to those obtained by NA49 and PHENIX collaborations, we keep most peripheral 
data in these two collisions out of the range of both the fits provided by eqn.(5) and eqn.(6). 

Obviously, the factors a and (3 provide the slopes of the fits. As the magnitude of any of 
these factors takes value nearer to zero, the fit will exhibit better fiatness of the data with 
respect to Nn-part or Ng^part implying a superior compliance of scaling , i.e., a better scaling. 
A look at Table-TVl reveals that the data on < tt >, and P/P in Pb + Pb collisions show 
a better degree of conformity with scaling when normalized by half of Ng_part, as in all cases 
1/3 1 is less than \a\. But, for the case of deuteron production in Pb + Pb collision, the scenario 
is just the opposite. Besides, the integrated yields for charged pions produced in Au + Au 
interactions at both the RHIC energies [Tables IVIl - 1 Vllj do not exhibit scaling satisfactorily, 
when the data are normalized by half of Ng^part instead of Nn-part, as the values of in 
these quite small with respect to |/?|. However, the data on K"^ and P/P produced 

in the same collisions favour the Nq_part-scalmg over the A^^-pa^^rScaling, though the cases 
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of P/P at ^/snn = 200 GeV are not as prominent as in the other cases. In fine, our net 
finding from the present study is the pions do not agree, while the kaons and protons do. 

The Fig. 6 (a) deserves some special attention and comments. The plots on integrated 
yields for charged hadrons produced in both nucleus-nucleus(y4 + A) and P -\- P collision 
defiect from each other when data on both are separately normalized in terms of number of 
participating nucleons(A^„_par.t)- The data on integrated yields in central nucleus-nucleus col- 
lisions normalized by the number of participant nucleon-pair have been used here from Fig. 3 
of Ref . [isjl . To our purpose,this has to be normalized in terms of the basic parton(quark)- 
constituents, denoted by < Ng^part >, for which the (charged)pseudo- rapidity density terms 
for various nucleus- nucleus collisions normalized by participant nucleon-pair, i.e., the factors 



5<N t> multiplied by R =< Nn-part > / < Nq^part > with the values of 

them as given in Table- IVTTTl And, in this conversion we have utilized those particular values 



of < Nn-part > wcrc uscd in Ref. jisl . so that we can normalize the exact values of 
by < Nq-part > /2. In calculating Ng-part for P + P collisions, we use equation(2) and the 
values of Ng-part for P -\- P collisions at three different cm. energies are given in Table- HXl 
function of relative probability. 
And when partonic considerations are used in normalization, the data on both P + P and 
A -|- A collisions come to an agreeable state. The < Ng-part >-values that are to be used for 
normalization of the most central P + P collisions, i.e.,0 — 5% central collisions are depicted 
in Table-El In order to check the nature of agreement between the data on A + A and P-\-P 
collisions, we try to obtain a fit by taking into account both sets of data, normalized by 
< Ng-part > /2; and the desired fit is to be described here by the following expression: 

1 dN,h _ _ (0.010 ±0.003) + (0.27 ±0.01) ln(^^) (7) 



0.5 < Ng_part > dr] GeV 
with /ndf = 3.632/12. The goodness of the fit is also shown in Fig. 6(b) with the help 
of the expression fit — data as a function of ^/Jnn which reveals that both sets of data are in 
good agreement with respect to the fit. Hence, both y4± A and P±P data exhibit a common 
yi^-dependence when the data on the former are normalized by the number of part.cipant 
parton-pairs. Very recently, Sarkisyan and Shakarov[I^ adopted a somewhat similar type 
approach in dealing with data on pseudorapidity density per participant-pair in cases of 
hadron-hadron, nucleus-nucleus collisions and also in the cases of e~^e~ interactions. The 
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study thus reinforced the idea of similarities of the bulk observables in the widest possible 
ranges of particle-interactions. 

IV. CONCLUDING REMARKS 

Let us now summarize our observations made here: (i) The secondaries, excluding 
deuterons, produced in Pb + Pb collision at CERN-SPS behave modestly well vis-a-vis 
Nq-part scaliug. The behaviour is more consistent towards the highest values of centrality, 
with centrality maximum at 0-5% (represented by solid boxes in the graph), (ii) Similar 
statements could be made about Au + Au interactions at both energies (Fig. 4 and Fig. 5) 
except the cases for pions. (iii) It is quite noticeable that there are modest degree of di- 
vergences in the A^„_par-rvalues between the calculations done by us at different centralities 
and those by others, as indicated in the text (in the second and third columns of the various 
Tables presented in this work). For the present, these discrepancies could be attributed only 
to our much simpler programming than what is resorted to by the big groups like NA49, 
WA98, PHENIX etc. The plain fact is we are now simply unable to take up such rigorous 
studies as might be desired or advisable due to various reasons beyond our capacity and 
control. But, we recognize the urgency and importance of such studies in order to arrive at 
a decision about the merit of Nq_part scaling and of the viewpoints expressed by Eremin and 
VoloshinQ. (iv) The Fig. 6 highlights on how the data on both P + P and A -|- A collisions 
negotiate better the idea of partonic-participant scaling, when the constituent-participant 
number for both P + P and A -|- A collisions are taken into account. Thus, the present 
study does essentially provide modest support to the viewpoints expressed by Eremin and 
Voloshin in their work of Ref . Q| . 
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TABLE I: 


Values of 


inelasitc 


nucleon-nucleon 


cross-section 


at different en- 


ergies. The second 


column 


provides the 


values obtained 


by subtracting 


eqn.(4) from 


eqn.(3). 


i lie ICXO L 


r'ol 1 1 1)1 11 (TTXTf^c; 
Uvjl LLiiili til V Co 


LliC VdlLiCo do 


L-lLt-U. Ill IVCi. tJ . 


/S AT AT(GeV) 






<Xf(mb) 
(Present Work) 




<^(mb) 
(From Ref.[5j) 


53 






35 




- 


56 






35 




37 


130 
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41 


200 






42 




42 


540 












630 






49 






900 






51 






1800 






56 








10 100 1000 10000 100000 

Vsnn (GeV) 

FIG. 1: Plots of total and elastic cross-section of PP and PP interactions 
as a function of cm. energies. Data points are taken from Ref. j^l. The 
solid curve depicts the fit for total cross-sections on the basis of eqn.(3) 
while the dashed one is for elastic cross-section on the basis of eqn.(4). 
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TABLE II: Values of < N^-part > and < Nq_part > for Pb + Pb collisions at E^ab = 
GeV. 



Centrality 




^ i\T PresentW ork ^ 
n~part ^ 


^ AT PresentW ork \ 


0-5% 


362 


368 


856 


5 - 14% 


305 


290 


644 


14 — 16/0 






44d 


ZA — 62/0 


1 on 

189 


14U 


278 


62 — 47% 


i3U 


no 


17i 


47 — 99%) 


11 


32 


53 


Centrality 


^ -'■^n—part 


^ 7y Presently orfc ^ 


^ q—part ^ 


0-1% 


380 ±1 


385 


901 


1 - 6.8% 


346 ± 1 


338 


770 


6.8 - 13% 


290 ±2 


273 


600 


13 - 25.3% 


224 ±1 


196 


411 


25.3 - 48.8% 


132 ±3 


99 


189 


48.8 - 67% 


63 ±2 


35 


56 


67 - 82.8% 


30 ±2 


10 


14 


82.8 - 100% 


12 ±2 


3 


4 
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TABLE III: Values of < Nn-part > and < Nq^part > for Au + Au collisions at ^/snn = 130 
GeV. 



Centrality 




^ AjPresentWork ^ 


^ AjPresentWi 
^ q~part 


0-5% 


347.7 ± 10 


351 


880 


5 - 10% 


293 ± 10 


297 


710.6 


5 - 15% 


271.3 ±8.4 


272 


645 


10 - 15% 


248 ±8 


242.5 


565.4 


15 - 20% 


211 ±7 


207 


469.4 


20 - 25% 


177 ±7 


173 


381 


15 - 30% 


180.2 ± 6.6 


174.5 


387.3 


25 - 30% 


146 ±6 


139 


298 


30 - 35% 


122 ±5 


120 


249 


35 - 40% 


99 ±5 


96 


192 


40 - 45% 


82 ±5 


81.5 


158.7 


45 - 50% 


68 ±4 


62 


114.6 


30 - 60% 


78.5 ±4.6 


77 


150.1 


60 - 92% 


14.3 ±3.3 


11.8 


17.8 
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TABLE IV: Values of < Nn-part > and < Nq^part > for Au + Au collisions at ^Jsnn = 200 



GeV. 



Centrality 






< 


iSjPresentWork ^ 
n — pavt 


^ AjPresentWork ^ 
^ q—part ^ 


0-5% 




351.4 ±2.9 




351.9 




880.1 


- 10% 




325.2 ± 3.3 




332.5 




819.7 


5 - 10% 




299.0 ± 3.8 




292.2 




094.0 


10 - 15% 




253.9 ±4.3 




242.4 




obU.y 


10 - 20% 




234.6 ±4.7 




226.2 




51 (A 


15 - 20% 




215.3 ±5.3 




208.2 




Add 

469 


20 - 30% 




166.6 ±5.4 




157 




339 


30 - 40% 




114.2 ±4.4 




108.4 




01 o 

zio 


40 - 50% 




74.4 ±3.8 




71.5 




1 Q /I Q 

134.3 


50 - 60% 




45.5 ±3.3 




40.4 




68 8 


60 - 70% 




25.7 ±3.8 




23 




35.6 


60 - 80% 




19.5 ±3.3 




16.6 




25 


fiO — 92% 

\J\J \J ^ / \j 




14.5 ±2.5 




11.8 




17.3 


70 - 80% 




13.4 ±3.0 




9.9 




1 3 6 


70 - 92% 




9.5 ± 1.9 




6.8 




Q 1 


80 - 92% 




6.3 ± 1.2 




3.9 
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TABLE V: 


Values 


of a, b, a and f3 for production of various 


secondaries : 


in Ph + Pb 


collisions at Ei^ah 


= 158A GeV[ndf= 


--No. of 


data - No. of parameters 


in the fit]. 


Identified 


a 


a 


X^/ndf 


b 


P 




Secondary 














< vr > 


0.74 ±0.06 0.13 ±0.01 


0.320/3 


0.82 ±0.08 


-(0.02 ±0.01) 0.442/3 




0.023 ±0.005 0.42 ±0.04 


3.436/3 


0.026 ± 0.006 


0.22 ±0.04 


3.700/3 


K- 


0.028 ±0.008 0.37 ±0.05 


11.024/3 


0.03 ±0.01 


0.18 ±0.05 


12.514/3 


P 


0.06 ±0.01 0.17 ±0.04 


4.492/3 


0.06 ±0.01 


0.03 ±0.01 


4.054/3 


P 


0.006 ± 0.002 0.24 ± 0.06 


4.616/3 


0.007 ± 0.002 


0.07 ±0.04 


4.313/3 


D 


0.0015 ± 0.0004 0.05 ±0.02 


0.635/3 


0.0017 ±0.0004 


-(0.11 ±0.04) 0.656/3 
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TABLE VI: Values of a, b, a and (3 for production of various secondaries in Au + Au 



collisions at ^/sNN = 130 GeV[ndf=No. of data - No. of parameters in the fit]. 



Identified 
Secondary 


CL 


Ci 


/ ndf 


b 


B 


/ ndf 




1.6 lb 0.1 


1 \ r\f\-\ r" 1 r\ r\f\r\'~) 

0.0015 ± 0.0003 


0.250/3 


1.6 ± 0.1 


— (0.13 ± 0.01) 


0.221/3 




1.3 ±0.1 


0.02 ±0.01 


0.190/3 


1.32 ±0.07 


-(0.12 ±0.01) 


0.125/3 


K+ 


0.11 ±0.01 


0.15 ±0.01 


0.073/3 


0.12 ±0.01 


-(0.012 ±0.002) 


0.083/3 


K- 


0.11 ±0.03 


0.11 ±0.06 


1.175/3 


0.11 ±0.03 


-(0.03 ±0.01) 


1.081/3 


P 


0.09 ± 0.01 


0.10 ± 0.02 


0.191/3 


0.09 ± 0.01 


-(0.05 ±0.01) 


0.151/3 


P 


0.06 ± 0.01 


0.10 ±0.02 


0.148/3 


0.06 ±0.01 


-(0.04 ±0.01) 


0.115/3 



TABLE VIL Values of a, 6, a and /3 for production of various secondaries in Au + Au 



collisions at ^/snn — 200 GeV[ndf=No. of data - No. of parameters in the fit]. 



Identified 
Secondary 


a 


a 


X^/ndf 


b 


/? 


X^/ndf 


7r+ 


1.7 ±0.1 


-(0.008 ±0.004) 


4.455/8 


1.6 ±0.1 


-(0.12 ±0.01) 


0.769/8 


7r~ 


1.4 ±0.1 


0.03 ± 0.01 


1.071/8 


1.5 ±0.1 


-(0.12 ±0.01) 


1.015/8 


K+ 


0.130 ± 0.004 


0.13 ±0.01 


0.583/8 


0.14 ±0.01 


-(0.03 ±0.01) 


0.591/8 


K- 


0.12 ±0.01 


0.14 ±0.01 


1.657/8 


0.13 ±0.01 


-(0.02 ±0.01) 


1.490/8 


P 


0.063 ± 0.002 


0.09 ±0.01 


0.264/8 


0.069 ± 0.003 


-(0.07 ±0.01) 


0.274/8 


P 


0.048 ± 0.002 


0.08 ± 0.01 


0.290/8 


0.053 ± 0.002 


-(0.07 ±0.01) 


0.347/8 
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TABLE VIII: Values of < Nn-part > and < Nq_part > used in 
Fig. 6 (a) to normalize ^^^-data for various nucleus-nucleus collisions 



Collision type 


y/SNN (GeV) 


Centrality 


TV T 7\ T ^ 

< ^^n~part > < ^^q-part > 


All + Au 


ALrb 


— 570 


3431151 786 


Ph + Ph 


8.7 


0-5% 


349 [16j 856 


Pb + Pb 


17.2 


0-5% 


362 [16] 856 


Au + Au 


56 


0-6% 


330 UJ] 823 


Au + Au 


130 


0-6% 


343 [12 871 


Au + Au 


200 


0-6% 


344[13] 871 


TART P lY- 
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FIG. 2: Comparison of average number of participant nucleons obtained in 
the present work and by different experimental collaboration's^ iii_ various li eav\ 
ion collisions. The experimental data are taken from Ref. 
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FIG. 3: Plots of integrated yields normalized by half of the number of participant- 
nucleons[(a),(b)] or constituent partons[(c),(d)] as a function of centralities for production 
of the various identified secondaries in Pb + Pb collisions at E^ab = GeVjsl, |l8|. The 

solid curves represent the fits obtained on the basis of eqn.(5)[(a),(b)] and eqn.(6)[(c),(d)]. 
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FIG. 4: Plots of integrated yields normalized by half of the number of participant- 
nucleons[(a),(b)] or constituent partons[(c),(d)] as a function of centralities for production 
of the various identified secondaries in Au + Au collisions at ^/sNN = 130 Gevjl^. The 
solid curves represent the fits obtained on the basis of eqn.(5)[(a),(b)] and eqn.(6)[(c),(d)]. 
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FIG. 5: Plots of integrated yields normalized by half of the number of participant- 
nucleons[(a),(b)] or constituent partons[(c),(d)] as a function of centralities for production 
of the various identified secondaries in Au + Au collisions at ^/s^n = 200 GevQ- The 
solid curves represent the fits obtained on the basis of eqn.(5)[(a),(b)] and eqn.(6)[(c),(d)]. 
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FIG. 6: (a) Energy dependences of integrated yields for charged hadrons produced 
in different nucleus-nucleus(74 + A) collisions at AGS, SPS and RHIC energies and 
for the same produced in P + P collisions at ISR energies(Fig.3 of Ref.Q]). The 
open boxes and open circles provide the data when normalized by half of the aver- 
age number of participant nucleons for nucleus-nucleus and PP collisions respectively. 
The solid boxes and solid circles show the same result for A + A and P + P colli- 
sions when normalized by half of the number of the participant-partons. The solid 
line provide a fit[eqn.(7)] for the nucleus- nucleus and proton-antiproton data when nor- 
malized by the half of the corresponding average number of constituent partons. (b) 
The plot of closeness of the different sejg of data with respect to the aforesaid fit. 
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